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The interaction between the cosmon and neutrinos may solve the "why now problem" for dark 
energy cosmologies. Within growing neutrino quintessence it leads to the formation of nonlinear 
neutrino lumps. For a test of such models by the integrated Sachs-Wolfe effect for the cosmic 
microwave background (CMB) we estimate the size and time evolution of the gravitational potential 
induced by these lumps. A population of lumps with size of 100 Mpc or more could lead to observable 
effects on the CMB anisotropies for low angular momenta. The linear approximation is found to be 
invalid for the relevant length scales. Quantitative estimates depend strongly on the details of the 
transition between the linear and nonlinear regimes. In particular, important backreaction effects 
arise from the nonlinearities of the cosmon interactions. At the present stage the uncertainties of 
the estimate make it difficult to constrain the parameter space of growing neutrino models. We 
explicitly discuss scenarios and models that are compatible with the CMB observations. 



o 

u 

43 

Oh- 
i 

o 



(N 
ON 

o 
o 



I. INTRODUCTION 

Quintessence models [2, Q with a growing neutrino 
mass can explain why dark energy dominates in the 
present cosmological epoch 0, A crucial ingredient 
is the coupling of the cosmon to neutrinos which results 
in an interaction between dark energy and matter [^, [|| . 
Indeed, the cosmon- neutrino coupling f3 can be large (as 
compared to gravitational strength) [7H9( , such that even 
the small fraction of cosmic energy density in neutrinos 
can have a substantial influence on the time evolution 
of the quintessence field or cosmon. This effect becomes 
important only once the neutrino mass plays a role. The 
transition from relativistic to nonrelativistic neutrinos at 
redshift znr ~ 5 — 10 therefore constitutes a "cosmologi- 
cal trigger event" which stops the further time evolution 
of the cosmon. For z < znr the time evolution of the cos- 
mic scale factor behaves very similar to an effective cos- 
mological constant which becomes computable in terms 
of the present day neutrino mass [H 0] ■ 

For realistic cosmologies the cosmon-neutrino coupling 
is larger than one, corresponding to an enhancement 
~ 1 of the cosmon-mediated attraction between non- 
relativistic neutrinos as compared to their gravitational 
attraction. In turn, this substantially accelerates the 
growth of neutrino structures for z < zjvjh, " t ne char- 
acteristic timescale is shortene d by a factor f3~ 2 . Large 
scale stable neutrino lumps 13, U| can form after a red- 
shift z n i Ri 1 — 2 at which the neutrino fluctuations be- 
come nonlinear [l2|, EH . The gravitational potential in- 
duced by neutrino lumps can have important effects on 
the propagation of photons and therefore influence the 
CMB-anisotropies via the integrated Sachs- Wolfe (ISW) 
effect. An extrapolation of the linear fluctuation anal- 
ysis in presence of the strong neutrino-cosmon coupling 
would indeed rule out most such models 14, 15j |. 



However, the linear analysis becomes invalid for the 
relevant length or momentum scales. An extrapolation 
of the linear evolution violates rapidly even the most ex- 
treme bound on the neutrino induced gravitational po- 
tential which would result if all neutrinos of the visible 
universe are concentrated in one spot The use of 

such an extrapolation for validity tests and parameter 
estimates of specific growing neutrino models [15| can 
yield misleading results. The challenge to be solved is 
an estimate of the neutrino-lump induced cosmological 
gravitational potential on large scales. 

The size of this potential and its time history deter- 
mines the size of the ISW effect which may be observed 
by a modification of the CMB-spectrum at low angular 
momentum I. It also influences the correlation between 
matter and photon fluctuations for which the present ob- 
servations yield values larger than expected in the cos- 
molog ical ACDM model [B-Eoj]. Finally, a rather sud- 
den increase of the neutrino-induced gravitational po- 
tential could reconcile [2l[ the presently observed large 
bulk flows [22j with observational bounds on the mat- 
ter fluctuations at similar scales [23|. If the attractive 
force between fluctuations is only mediated by gravity 
rather tight bounds relate the growth rate and the pecu- 
liar velocities to the density contrast. This bound gets 
weakened if the time scale is shortened, as characteristic 
for the nongravitational cosmon-mediated force. 

The difficulty of a reliable estimate of the neutrino- 
induced gravitational potential comes from the simple 
observation that a gas of neutrino lumps behaves dif- 
ferently from a fluid of individual neutrinos. This con- 
cerns both gravitational and cosmon-mediated effects. 
On small enough length scales, where neutrino lumps 
form first, the nonlinear part of the evolution is domi- 
nated by the formation of individual lumps. We have 
computed the time history of individual lumps by solv- 
ing the relevant nonlinear Navier-Stokcs equations in the 
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Newtonian limit 13|. The local gravitational potential 
on a length scale characteristic for the lump increases 
faster than in the linear approximation during an "infall 
period". Subsequently, it settles to an almost constant 
value. With reasonable assumptions about the distribu- 
tion of the lumps these results can be used for an estimate 
of the neutrino induced cosmological potential $ v (k). 
The neutrino fluctuations in this "nonlinear range" of 
rather large k have no substantial effect on the CMB. 

On the other hand, for length scales sufficiently large 
compared to the horizon, linear perturbation theory ap- 
plies. Again, the neutrino fluctuations in the "linear 
regime" have only minor effects on the CMB. Visible 
effects can arise from neutrino induced fluctuations of 
the gravitational potential $ v (k) in the intermediate 
range of scales, typically only moderately inside the hori- 
zon, which interpolates between the linear and nonlinear 
ranges. A potential $„ (k) exceeding substantially a value 
around 10 -5 would have a large impact on the CMB and 
rule out the corresponding model. 

The gravitational side of the breakdown of linear per- 
turbation theory for neutrino fluctuations resembles in 
some aspects the situation for cold dark matter fluctua- 
tions, but also shows important differences. The latter 
is well understood. Cold dark matter fluctuations grow 
nonlinear first on scales much smaller than the horizon. 
For example, galaxies are formed long before the emer- 
gence of superclusters, happening in the present cosmo- 
logical epoch. The complicated nonlinear formation pro- 
cess of a galaxy has only little impact on its gravitational 
potential at a distance I if the infall of additional mat- 
ter into a volume Z 3 around the galaxy can be neglected. 
Once the infall of new matter stops, the gravitational po- 
tential of a single galaxy for a characteristic scale I does 
not change in time anymore. The cosmological gravita- 
tional potential (more precisely the root of the squared 
fluctuations of the potential in Fourier space) at a co- 
moving wave number k = na/l changes then dominantly 
due to the expansion of the universe and decreases due 
to dilution. 

In principle, the formation of galaxies could influence 
the growth of structure on length scales much larger than 
galactic scales. Beyond linear perturbation theory the 
fluctuations with different k are coupled such that a for- 
mal expansion in lowest order of the fluctuations breaks 
down once the galactic-scale fluctuations grow nonlinear. 
In principle, a gas of galaxies needs not to be the same 
as a fluid of dark matter particles. This extends to larger 
nonlinear cosmic structures like clusters or local voids. 
For the case of gravity, however, the question if structure 
formation can influence the cosmological evolution has 
been answered to the negative 24j . This applies similarly 
to the evolution of fluctuations for k much smaller than 
the characteristic k n \ where fluctuations have grown non- 
linear. The central point of this argument states that the 
gravitational field sufficiently outside a mass concentra- 



tion or energy density inhomogeneity is quite insensitive 
to the details of the mass distribution or local relative 
velocities. The gravitational field equation is linear in 
the energy momentum tensor such that nonlinear "back- 
reaction" effects can only arise due to self-interactions of 
the gravitational field. For realistic cosmological struc- 
tures a typical estimate of the gravitational backreac- 
tion yields relative effects only of the order 10 -5 [2^ - [26j | . 
In short, linear perturbation theory for fluctuations at a 
scale k is not affected by nonlinearities at larger k' 3> k, 
provided that the gravitational potential $(fc') induced 
by these "local lumps" remains small compared to one. 
(Challenges of this simple picture by s pec ulations about 
large gravitational backreaction effects [27] appear not to 
be compatible with a realistic distribution of structures 
in the universe). The decoupling of the small distance 
nonlinearities from the linear evolution of large distance 
inhomogeneities is the basis for the predictions of the 
CMB-anisotropies in the ACDM model. 

For neutrino fluctuations the gravitational aspects of 
the difference between a gas of neutrino lumps and a 
fluid of neutrino particles differ from cold dark matter in 
two important aspects. The first concerns the difference 
in the relevant length scales. For cold dark matter the 
scales of substantial nonlinearities do not exceed much 
the scales of clusters and the nonlinear Rees-Sciama ef- 
fect affects the CMB only for angular momenta / > 1500 
[28| . Due to the additional cosmon-mediated attrac- 
tive force, the neutrino fluctuations become nonlinear 
for much smaller k as compared to dark matter. Non- 
linear lumps with size exceeding considerably 100 Mpc 
may form. This implies that the breakdown of linear per- 
turbation theory affects the CMB-anisotropies for much 
smaller angular momenta I. For the models that we dis- 
cuss in this paper the linear approximation for the CMB 
breaks down for I < 200; for k > 10 Mpc -1 the neu- 
trino fluctuations remain linear due to the effect of free 
streaming. 

The second difference is related to the dependence of 
the neutrino mass on the value of the cosmon field. Its 
value inside the lumps may differ substantially from the 
average cosmological value which would be computed for 
an almost homogeneous neutrino fluid. Since the mass is 
the source of the gravitational potential at all scales, a 
concentration of a substantial fraction of neutrinos within 
lumps where their mass differs from the cosmological 
value could result in important modifications for the evo- 
lution of the gravitational potential on large length scales 
due to nonlinearities on small scales. This possibility of a 
large backreaction effect differs from pure gravity, where 
the particle masses are constant. 

The main differences between a neutrino fluid and a 
gas of neutrino lumps concern the cosmon field rather 
than the gravitational aspects. For growing neutrino 
quintessence the dominant attractive force between non- 
relativistic neutrinos is mediated by the cosmon. In a 
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Newtonian approximation the role of the gravitational 
potential is taken by /3 54>, with 5<p the local fluctuation 
of the cosmon or quintessence field. While the gravita- 
tional potential remains much smaller than one even for 
highly nonlinear neutrino lumps, the quantity /3 8<f> can 
reach values of order one, signalling a breakdown of linear 
theory in the cosmon sector. Indeed, in a Newtonian ap- 
proximation (3 8(j) exceeds the gravitational potential by 
a factor 2f3 2 which reaches values 10 4 - 10 6 in our mod- 
els. A local value (3 5<fi = 1 means that the local mass of 
the neutrinos is smaller by a factor e _1 as compared to 
the cosmological neutrino mass. One expects substantial 
modifications of linear theory once (5 5<p reaches values of 
order unity. 

In the presence of a cosmon field the impact of small 
scale structure on the behavior of fluctuations at larger 
scales (backreaction) is generically expected to be much 
larger than in the case of gravity 24, [29[. The rea- 
son are the nonlinearities in the cosmon field equations. 
They concern both the role of the cosmon potential and 
cosmon-neutrino coupling, in particular if the effective 
interaction strength /3(0) depends on the value of the cos- 
mon field. As mentioned above, these nonlinearities may 
cause a difference between the average neutrino mass in- 
side the neutrino lumps as compared to the value which 
would be obtained for a roughly homogeneous fluid of 
neutrinos. Also the effective cosmon-neutrino coupling 
p = — dh\m v / dej) may be reduced for a gas of neutrino 
lumps as compared to a homogeneous neutrino fluid. 

In this paper we account for the effect of the nonlin- 
earities by simple physically motivated bounds on the 
growth of fluctuations. This will allow us to give a first 
estimate of the effects of neutrino lumps on the CMB 
anisotropics. Our findings clearly indicate that linear 
perturbation theory is not appropriate. On the other 
hand, the uncertainties of the quantitative treatment of 
the nonlinear effects remain large, such that reliable es- 
timates of the CMB spectrum are not yet available for 
growing neutrino quintessence. The latter would require 
a better quantitative estimate of the cosmological grav- 
itational potential induced by neutrino lumps, perhaps 
by numerical N-body simulations or nonperturbative an- 
alytic methods. 



II. GROWING NEUTRINO QUINTESSENCE 

The time evolution of growing neutrino models follows 
a set of equations .12] which combines the evolution of the 
homogeneous and isotropic background with the growth 
of perturbations. We recall for convenience the essential 
ingredients characterizing these models. At the back- 
ground level, the expansion of the universe obeys the 



Friedmann equation: 
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(i) 



where primes denote derivatives with respect to confor- 
mal time r. The sum is taken over all components a of 
the energy density in the universe, including dark mat- 
ter, quintessence, neutrinos, baryons and radiation. The 
time evolution of the energy density p a for each species 
involves the equation of state w a = p a /p a - 

A crucial ingredient in this model is the dependence 
of the average neutrino mass on the cosmon field 0, as 
encoded in the dimensionless cosmon-neutrino coupling 



d\nm l , 



(2) 



For increasing <p and constant /3 < the neutrino mass 
increases with time 



m v e 



-P<S> 



(3) 



where m„ is a constant. In general, /3 can be a func- 
tion of cf>, as proposed in Q within a particle physics 
model, modifying equation ([3]) but leading to similar ef- 
fects. The cosmon field <f> is normalized in units of the 
reduced Planck mass M = {SttGn)- 1 / 2 . For a given cos- 
mological model with a given time dependence of <f>, one 
can determine the time dependence of the neutrino mass 
m u {t). For three degenerate neutrinos the present aver- 
age value of the neutrino mass m v (to) can be related to 
the energy fraction in neutrinos (h « 0.72) 



a, (t ) 



MeVh 2 



(4) 



The cosmon exchange mediates an additional attrac- 
tive force between neutrinos of strength f3 2 . The case 
j3 1 corresponds to a strength comparable to gravity. 
In this paper we investigate two models, one with con- 
stant (3 = —275 and m v (to) = 0.47 eV, and the other 
with varying /3(0) as proposed in In this "growing 
coupling model" an approach towards zero of the large 
mass scale in a see-saw type mechanism for the neutrino 
mass is reflected by a strong increase of the light neutrino 
masses that we take all equal for simplicity, 



m v <p t 



1 



, (5) 



with constants m v and 4>t- In a coupled model this be- 
havior is only approximate. The pole at <f> = <j>t is not 
reached in practice. We choose a present neutrino mass 
m u (t ) = 0.46 eV and /3(t ) = -150. 

The dynamics of the cosmon can be inferred from the 
Klein Gordon equation, now including an extra source 
due to the neutrino coupling, 
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+ = a 2 13{p v - i Pu ) , (6) 



with p v and p v — w vPv the energy density and pressure 
of the neutrinos. We choose an exponential potential 
ASM, 0: 



V{(j>) = M 2 U((j)) = M 4 



(7) 



where the constant a is one of the free parameters of our 
model and determines the amount of early dark energy. 
Current bounds constrain it to be of the order a ~ 10 
or bigger [32| and here we use a — 10. The homoge- 
neous energy density and pressure of the scalar field <fi 
are defined in the usual way as 



2a? 



+v(4>) , P<t> = 



2a 2 



-V{(j>) , w<j, 



(8) 



We choose the value of 7 = — (3 /a such that we obtain 
the correct present dark energy fraction according to the 
relation [H, 



fift(to) = 



16eV 



7(*o) 



(9) 



Finally, we can express the conservation equations for 
dark energy and growing neutrinos as follows 0, Q : 

^ = -3H(l + w W+/V(l-3uvK , (10) 

P ' v = -m{i + w v ) Pv - p<t>'{i-zw v )p v . (11) 

The sum of the energy momentum tensors for neutri- 
nos and dark energy is conserved, but not the separate 
parts. We neglect a possible cosmon coupling to Cold 
Dark Matter (CDM), so that p' c = -3Hp c - 

Given the potential ([7]), the evolution equations for the 
different species can be numerically integrated, providing 
the background evolution shown in Fig|T] for the model 
with constant j3. The initial pattern is a typical early 
dark energy model, since neutrinos are still relativistic 
and almost massless, withp„ = so that the coupling 

term in eqs.(j6|), (fT0|) . (fTTj) vanishes. Dark energy is still 
subdominant and falls into the attractor provided by the 
exponential potential (see [j||6|,[33| for details). For z > 5 
it tracks the dominant background component with an 
early dark energy fraction fl^ = n/a 2 and n — 3(4) for 
the matter (radiation) dominated era. 

As the mass of the neutrinos increases with time, the 
coupling term ~ j3p v in the evolution equation for the 
cosmon ([6]) starts to play a significant role as soon as neu- 
trinos become nonrelativistic. This cosmological "trigger 
event" kicks <j> out of the attractor. For z < 5 the value of 
the cosmon stays almost constant, and the frozen cosmon 
potential mimics a cosmological constant. In Fig|T]this 
is visible in the modified behavior of p u and for z < 5. 



Subsequently, small decaying oscillations characterize the 
<j}—v coupled fluid and also p v reaches an almost constant 
value. For constant j3 the increase of m u ~ a 3 cancels the 
dilution of the neutrino number n v ~ a -3 . This feature 
is altered for the growing coupling models (FigJ5|). The 
values of the energy densities today are in agreement with 
observations, once the precise crossing time for the end 
of the scaling solution has been fixed by an appropriate 
choice of parameters. At present the neutrinos are still 
subdominant with respect to CDM, though in the future 
they will take the lead. 
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FIG. 1. Background evolution for fixed cosmon-neutrino cou- 
pling. Energy densities of cold dark matter (solid, red), neu- 
trinos (long-dashed, green), dark energy (dot-dashed, blue) 
and photons (short dashed, black) are plotted vs redshift, for 
P = -275, a = 10, m v (t ) = 0.48 eV. 



III. LINEAR AND NONLINEAR EVOLUTION 
OF NEUTRINO FLUCTUATIONS 

Once neutrinos become nonrelativistic (znr ~ 5), 
they start feeling the attractive force stronger than grav- 
ity mediated by the cosmon field. Neutrino perturba- 
tions rapidly grow and become nonlinear at a redshift 
z u \ ^1 — 2 [121 ] . Subsequently they might form stable 
lumps, whose solutions have been described in 
This scenario poses a challenge for our understanding of 
the spectrum of CMB anisotropies in the range of angu- 
lar momenta I ;$ 200. In dependence on the size of the 
gravitational potential generated by the neutrino lumps 
one may find a strong influence on the ISW effect: pho- 
tons falling into a potential well may come out with a 
different energy if the size of the gravitational potential 
has changed in the meantime. Details depend on the 
size and time evolution of the gravitational potential at 
various scales k. 
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FIG. 2. Background evolution for varying cosmon-neutrino 
coupling. Energy densities of cold dark matter (solid, red), 
neutrinos (long-dashed, green), dark energy (dot-dashed, 
blue) and photons (short dashed, black) are plotted vs red- 
shift, forgrowing coupling, a — 10, m„(to) = 0.46 eV. 



In the relevant range of k and z the linear fluctua- 
tion equations for the neutrinos and the cosmon break 
down. Any realistic estimate for the size of the modi- 
fications of the CMB-spectra needs an estimate of the 
neutrino induced gravitational potential in the nonlin- 
ear regime. An extrapolation of the linear fluctuation 
equations .15] would give a large impact on the CMB- 
anisotropies which would rule out the models discussed 
in this paper. We will see, however, that this effect is 
overestimated by many orders of magnitude. In this pa- 
per we attempt a first estimate of the neutrino induced 
gravitational potential, based on extrapolations relying 
on the linear equations for small k < 10~ 4 h/Mpc and 
nonlinear results for k > 2 • 10 -2 h/Mpc. 

A summary of our findings is shown in Fig|3l where we 
display various approximations to the neutrino-induced 
gravitational potential $ v (k) for the present cosmological 
time (z — 0). Here $„(fc) is related in the Newtonian 
approximation to the neutrino density contrast A u 



p v {z)A v {k,z) 
2k 2 {l + z) 2 



(12) 



where the factor (l+z)~ 2 — a 2 (z) arises from the relation 
between comoving wave number k and physical momen- 
tum q — k/a. The neutrino density contrast in Fourier 
space A„(fc) is obtained from the two point function or 
power spectrum P v (k) = V(\S„(k)\ 2 ) as 



A 2 (fc) 



k z P„{k) 



S v (k) = V- L I d 6 x8 v (x)e~ 1 ^ , S v (x) 



Pv 



(14) 

In this figure, we have plotted the linear estimate (rising 
curve on the left) together with simple non-linear sce- 
narios (falling curves on the right) and interpolations. 
Origin and details of these estimates will be given later 
in the text. 

Unfortunately, the effects on the CMB spectrum are 
dominated by the gravitational potential for scales near 
the horizon, k « 3 • 1(T 4 h/Mpc. This is the region 
where the interpolation suffers from large uncertainties. 
As an overall picture we may state that interesting mod- 
ifications of the CMB spectrum for I ;$ 60 occur if the 
neutrino induced gravitational potential reaches a value 
<I>„(fc) of a few times 10~ 6 . For much larger values of 
<!>„ the models are ruled out, whereas for much smaller 
values the effects on the CMB become insignificant. 

For small enough k and large enough z the linear fluc- 
tuation equations remain valid. For our system of cou- 
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(13) 



pled fluids they can be found in ref. [12| ■ The estimate of 
the scales and redshifts where the linear equations break 
down is not easy. There is an obvious breakdown when 
the neutrino overdensity A v reaches unity. In practice, 
the critical value A£ for the breakdown of the linear evo- 
lution may be considerably smaller. A second impor- 
tant quantity is the product /3 S4>, where 8<f> denotes the 
fluctuation of the cosmon field around its cosmological 
background value. For local values \/3 S4>\ s=s 1 the lo- 
cal neutrino masses differ strongly from the cosmological 
background value. The linear approximation assumes im- 
plicitly that the neutrino masses are given by the back- 
ground value. It breaks down if this does no longer hold 
for a substantial fraction of the neutrinos. This extends 
to the effective value of the cosmon neutrino coupling ft. 

The hardest problem in this context concerns the issue 
of "backreaction" of "small scale fluctuations" well in- 
side the horizon on "large scale fluctuations" close to the 
horizon. Evaluated for a given local neutrino lump the 
nonlinear effects arising for |/3<5</>| ^1 1 occur already for 
rather large z. Subsequently, the neutrino mass within 
a virialized lump may differ from the cosmological value 
outside the lump. As an effect, the fluctuations on larger 
length scales may "see" a modified effective neutrino 
mass and a modified ^-dependence of m v {_<t>). This back- 
reaction of fluctuations with larger k on the evolution of 
fluctuations with smaller k tells us that the bound of va- 
lidity for the linear fluctuations at a given scale k may 
not only depend on quantities evaluated at this scale. 

A local value j35(j>^ w 1 for a lump of size I corresponds 
to a much smaller value of the cosmological value f35<j)(k) 
with k = iva/l. The averaging over the horizon volume 
contained in the Fourier transform leads to a suppression 
factor 7 C which is given roughly by the total volume of 
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all lumps with size ~ I divided by the horizon volume 
6<j>(k) « j c 54>i(l = ira/k). In [13[ this suppression factor 
has been estimated as j c s» 10~ 3 . Thus linear pertur- 
bation theory could break down already for rather small 
values j3S<f) s» 10~ 3 . An explicit analysis of higher order 
terms in sect fVl shows that this is indeed the case. 

On the opposite side for large k and small z the lin- 
ear equations have clearly become meaningless. In order 
to understand the evolution of single neutrino lumps on 
scales well within the horizon one can solve numerically 
(HI the Navier- Stokes equations in the Newtonian ap- 
proximation. In an expanding universe with comoving 
coordinates x they read in position space 



< = 00') v„-(v„-V)v„ 

AScp = —/3 a 2 8 v p v , 



(15) 

(16) 
(17) 

(18) 



Here p v is the background neutrino energy density and 
8 V = Spv/pv is the relative neutrino density perturba- 
tion. The vector v„ is the velocity for neutrinos. More 
precisely, it describes the peculiar comoving velocities - 
it vanishes for neutrinos with constant comoving coordi- 
nates. The evolution of the velocities is driven by the 
gradients of the gravitational potential and cosmon field, 
with the usual Hubble damping and quadratic term aris- 
ing from particle number conservation. The velocity de- 
pendent term /30'v„ in eg. (1161) is not present in the stan- 
dard Navier Stokes equations. It accounts for momen- 
tum conservation, reflecting the fact that the neutrino 
mass changes in time as m'/m = —/3<f>'. Equation (|18l) is 
the Poisson equation for the gravitational potential that 
we have indexed as to clarify that it only comprises 
the neutrino contribution. Correspondingly, ea.(|17[) de- 
termines the local scalar potential 5(f). In this equation 
we have neglected the contribution from the cosmon po- 
tential U{(f>). For the growing coupling model cq. (fTT|) 
receives an additional contribution ~ df3/d(f> that we de- 
scribe in sec|Vl We further note that the gravitational 
potential in position space as defined in eq. (|18[) is not 
the Fourier conjugate of $„(fc, z) in eq. (|T2j) . due to the 
definition of A(fc) in eq. ([T^]). For a detailed discussion 
of the underlying equations and the numerical solution 
techniques we refer to [H. 

In ref. [13| it was found that neutrino lumps virialize 
at a redshift z V i r ps 1.2 for m„(£o) = 2.1 eV and a typical 
size Rf = 14 Mpc of the lump. Here we define the scale 
R f roughly as the scale where the gravitational potential 
reaches about one half of its value in the inner region of 
the lumps at the time of virialization. For single lumps it 
seems reasonable that the gravitational potential at the 
length scale Rf does not change substantially for z < 
z v i r , since the details of what happens in the core of the 



lump should not matter for gravity outside. 

The extrapolation from the gravitational potential of 
a single lump to the average cosmological gravitational 
potential at the corresponding scale k is difficult. Since 
the latter involves a Fourier transform over distributions 
of lumps of different sizes, an estimate of the neutrino 
induced gravitational potential $ u (k) for strongly non- 
linear scales has to assume some estimate of the number 
and distribution in size and mass of the virialized neu- 
trino lumps. At the present time this estimate involves 

Nevertheless, some lessons 
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important uncertainties 
on the scale dependence, z-dependence and absolute up- 
per bounds for <fr u (k) can be drawn. We will discuss this 
issue in more detail in the next section. 

The most uncertain part of our estimate of the fluctua- 
tion effects concerns the interpolation between the linear 
fluctuations and the "individual lump regime" . This is 
the regime where the dynamical behavior of a gas of neu- 
trino lumps matters and where backreaction effects could 
be important. This interpolating regime makes the dom- 
inant contribution to the ISW effect. 

After these general cautionary remarks we may now 
turn back to FigJ3] where we show the neutrino induced 
cosmological gravitational potential <fr v (k) at z — for 
the different regimes in fc, evaluated for the model with 
varying f3. The extrapolation of the linear perturbation 
equations (green long-dashed line) leads to a strong in- 
crease of $iy(fc) even for scales somewhat outside the hori- 
zon, where in absence of the cosmon mediated attraction 
between inhomogeneities one would have Q v ~ k~ 2 . As 
a consequence, the linearly extrapolated <fr v (k) hits the 
extreme bound (for which all neutrinos within the hori- 
zon are concentrated in one spot) for length scales close 
to the horizon. As k increases it soon reaches values that 
correspond to the potential near the Schwarzschild radius 
of a black hole. It is obvious that this is not meaningful. 
We have indicated by a solid dot on the dashed green 
line the value of k for which the linear approximation to 
A„(fc) reaches the value 1/4, with larger values of A w 
for larger k. In other words, k n i = 1.71 • 10~ 3 h Mpc" 1 
(corresponding to R n i = hir/k = 2.6 Gpc) marks the 
smallest length scale for which the final linear density 
perturbation obeys |A„(fc)| = 1/4. One sees that neu- 
trino fluctuations on scales only slightly smaller than the 
horizon, Rh ~ 3 Gpc/h, are already highly nonlinear. 

In the region k > 0.01 we indicate in Fig |3] the result of 
rather radical and oversimplified assumptions. We sup- 
pose that one quarter of all cosmic neutrinos is clumped 
in lumps of a given size. For each lump we take the 
density profile that has been found from the solution of 
the hydrodynamical nonlinear equations for individual 
lumps in [13|. We display two different sizes of lumps, 
Rf = 15 Mpc and 150 Mpc. We expect the true neutrino 
induced gravitational potential to be below these curves. 
The decrease of the curves for large k is due to the finite 
"effective width" of the neutrino distribution within the 
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lumps. Concentrating all neutrinos of the lump in the 
center yields the dotted straight lines in Figj3] The ex- 
treme bound for <Ev(fc) assumes that all neutrinos within 
our horizon are concentrated in one point. It is indicated 
by the black dashed line. Indeed, if all neutrinos within 
the horizon volume Vh or are concentrated in one point 
one finds 

(19) 

If all neutrinos arc distributed in TV equal pointlikc 
lumps, the gravitational potential (|19[) is multiplied by 
a factor N^ 1 / 2 , and a further factor f v arises if only a 
fraction f v of the neutrinos within the horizon are dis- 
tributed in such pointlike lumps. 

In the absence of nonlinear effects in the cosmon sec- 
tor a typical interpolation would follow the line of linear 
growth until it hits a critical value, for which we take the 
extrapolation of eq. fH?)) with A„(fc, z) = A cr jt- In FigJ3] 
we indicate the value A cr ;t = 1/4 by a dot. Even though 
one expects some smoothcning in the region where the 
two curves match, the resulting values of typically 
would exceed a value 10 -5 in the matching region. The 
ISW-contribution would in this case exceed the observa- 
tional bounds. 

For the interpolation shown in Fig(3] the increase of 
$ v {k) outside and close to the horizon is damped by back- 
reaction effects. This interpolation joins rather smoothly 
with the nonlinear regime for individual neutrino lumps 
for k > 2 ■ 10~ 2 h Mpc -1 . The large values of <& u (k) close 
to the dot are never reached. The upper one of the blue 
dotted curves leads to an observable ISW-effect, while 
for the lower one no substantial impact on the CMB is 
found. This demonstrates the sensitivity of a possible 
CMB-signal to the details of the interpolation. 



IV. GRAVITATIONAL POTENTIAL OF 
NEUTRINO LUMPS 

The local gravitational potential of a single lump can 
be directly accessed via a numerical integration of equa- 
tions (rT5|) - (IT8l) 13] . In Fig|4]we plot the time evolution of 
the gravitational potential of a single neutrino lump with 
final radius Rf — 14.3 Mpc for different growing neutrino 
models. Here Rf is the distance from the center where 
the gravitational potential reaches one half of the value in 
the core at virialization. We show the local potential (in 
position space) (Rf). In either model, the collapse of 
the lump and the growing neutrino mass lead to an in- 
crease of the gravitational potential. The 0- dependence 
of m„ is responsible for the observed small oscillations. 
Two features distinguish the particular models consid- 
ered: 



(i) Onset of growth: In the model with a large con- 
stant p — —275 the growth of the neutrino mass 
is very rapid. Due to the smaller mass m v (to) = 
0.47 eV, the neutrinos become nonrelativistic later 
(znr ~ 2) than in the models with constant /? = 
—52 and large neutrino mass m„(io) =2.1 eV, for 
which one finds znr ~ 5. For the growing coupling 
model one has zjvi? ~ 3. For smaller values of znr 
the growth of structure sets in later. For the red 
solid curve (/3 = —275) the neglection of pressure 
terms becomes justified only for z < 2. 

(ii) Time until virialization: On the other hand, the 
stronger coupling leads to a substantial increase of 
the growth rate in the model with /3 = —275. Hence 
the time it takes from approximate linearity to viri- 
alization is much shorter in this model (Az ~ 0.4) 
as compared to Az ~ 3 for the model with constant 
(3 = —52 and large neutrino mass m„(to) = 2.1 eV. 
This time is even longer for the varying coupling 
model. 

The combination of both effects leads to a rather similar 
redshift of virialization z v j r for the different models. For 
the two models under investigation we expect virialized 
structures around redshift z v ; r ~ 1.2. 

For an estimate of the cosmological impact of neutrino 
lumps one needs the average cosmological value of the 
gravitational potential at a given scale, expressed as the 
Fourier component $j,(fc). By distributing lumps of given 
sizes in a cosmological volume, we have obtained rough 
estimates of this quantity as a function of the comoving 
scale k. We refer to [l3| for a detailed description of the 
procedure. As a rough approximation to these non-linear 
gravitational potentials in the range 0.03 Mpc -1 < k < 
0.3 Mpc -1 after virialization we find 

$ u (k,z) = c{z)k- a , aw -1.6 . (20) 

For the rather crude nonlinear scenarios shown in Fig|3] 
we obtain c(z V i r ) ~ 1.2 x 10 -8 , with k in units of Mpc -1 . 

The curves in Fig J3] reflect the result of the rather rad- 
ical assumption that all lumps form without substan- 
tial disturbance from other lumps in a situation close 
to spherical infall. Furthermore, the cosmological dis- 
tribution of such lumps is strongly oversimplified by the 
assumption that a fraction /„ = 1/4 of the cosmic neu- 
trinos is found in nonlinear lumps, and that all lumps 
are of equal size. We expect several effects to reduce 
the corresponding value of the neutrino induced gravi- 
tational potential <& v (k). A first substantial reduction 
arises if scattering and merging of smaller lumps plays a 
dominant role for the formation of larger lumps, rather 
than the roughly spherical infall for the formation of sin- 
gle lumps. This holds, in particular, if one attempts to 
extrapolate the estimate of the nonlinear $„(fc) to scales 
k 6.3 • 10~ 3 h/Mpc where the Newtonian approxima- 
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FIG. 3. Neutrino induced gravitational potential as a func- 
tion of comoving wave number k for 2 = 0. We show re- 
sults for the growing coupling model with present average 
neutrino mass m v (to) = 0.46 eV. The red and light blue lines 
are inferred from the nonlinear evolution. The solid lines in- 
dicate simple scenarios where one quarter of all neutrinos are 
clumped in extended lumps with the same size of 15 (red) and 
150 Mpc (blue), respectively. The dashed lines are obtained 
from lumps concentrated in one point of corresponding mass. 
The behavior of <E>„ according to the linear approximation is 
shown by the long-dashed, green line. The dot indicates the 
value of k for which the linear extrapolation of A„(fc) reaches 
the value 1/4. The dot-dashed black line shows the extreme 
bound (|19[) . The interpolating dotted lines correspond to the 
bounded linear gravitational potential with rj cr i t = 0.001 and 
rjcrit = 0.01, and A crit = 0.1 and A cr i t = 1/4, respectively, 
see details in sec. (I VI t). 



tion breaks down. Furthermore, backreaction effects typ- 
ically slow down the growth of larger size neutrino struc- 
tures, again reducing the cosmological value |$„(fc)| as k 
decreases (cf. the interpolated curve in FigJ3]). 

Finally, the neutrino-induced gravitational potential 
after virialization grows ~ m v (z). This follows if we as- 
sume that the neutrino number distribution n v {r), as 
a function of physical distance r from the center of the 
lump, does not change any further after virialization. For 
the nonlinear estimates we have neglected possible back- 
reaction effects and use the cosmological neutrino mass 
m v {z) according to the background evolution J2]), ([3]). 
This typically leads to an increase of $„(fc) by a fac- 
tor (1 + Zvir) 3 ~ 10 after virialization. The true value 
of the gravitational potential of neutrino lumps should 
rather use the average value of the neutrino mass in the 
lumps which may be substantially smaller than the cos- 
mological value. In contrast, the lower two interpolated 
curves effectively stop the growth of m v within the lump 
once nonlinear effects become important. After virializa- 
tion the time evolution of <& v (k) arises then only from the 
rescaling between comoving and physical momenta and a 



FIG. 4. Local neutrino gravitational potential as a func- 
tion of redshift for a fixed final effective gravitational radius 
Rf = 14.3 Mpc of the lump. We show models with fixed 
/3 = -275 (solid, red) p = -52 (dotted, blue) with corre- 
sponding neutrino masses m v {to) = 0.47 or 2.11 eV, respec- 
tively. The long-dashed, green curve refers to the growing 
coupling model with m v (t ) — 0.46 eV. 



factor a -3 / 2 from cosmological averaging (cf. sect. (jVI|) ), 
This discussion reveals quite large potential errors in our 
nonlinear estimates which may rather be interpreted as 
approximate upper bounds. A reduction of <& v {k) by a 
factor 10 or more seems quite reasonable, especially for 
small fc. 

It is also worth to understand the evolution and break- 
down of linear perturbation theory in more detail. In the 
linear approximation the gravitational potential is given 
by the expression 



a 2 / H 



(21) 



The first term denotes the density contribution and 
should match the nonlinear density contribution after 
virialization, whereas the second term, the "velocity con- 
tribution" approximately vanishes after virialization. In 
Fig JS] we plot the linear evolution for 4> u (q), evaluated at 
the physical scale q = k/a — 6.3 ■ 10~ 3 h/Mpc. This 
should never exceed upper bounds obtained from the 
nonlinear evolution. In Figf5]we compare this linear ap- 
proximation with the crudely oversimplified assumption 
that only neutrinos in randomly distributed single lumps 
of final size Rf = ir/q contribute, and that a fraction 
fv = 0.25 of the total number of neutrinos is bound in 
such lumps. This results in the almost horizontal line in 
the left part of FigfSJ We may regard the corresponding 
green line as an upper bound for $„ at this scale. What is 
clear at this stage already is the huge overestimate of $„ 
which would result from the extrapolation of the linear 
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approximation to the density contribution to Q v (shown 
in Figj5] as the solid red line) . 



Linear <£ v 
Nonlinear <P., 




1+z 



FIG. 5. Breakdown of the linear approximation for the grav- 
itational potential for a physical wave number q = k/a — 
6.3 ■ 10" 3 h/Mpc. For the growing coupling model we show 
the redshift dependence of the linear approximation of <& v 
as the solid red line. It exceeds the rather extreme upper 
bound for from a nonlinear estimate for final lumps of ra- 
dius Rf = 15 Mpc which is given by the horizontal (wiggled) 
green line for redshifts after virialization. 



In FigJBl we show for the growing coupling model the 
linear approximation to <fr u (k) for various redshifts. This 
is confronted with randomly distributed nonlinear lumps 
of given size, represented as a power law approximation 
to the (scattered) nonlinear curves in FigJ3] We also in- 
dicate the extreme bound where all neutrinos are concen- 
trated in one point. For z = 3.4 the linear approximation 
remains valid for the entire range of k. For z = Zj, = 1.7 
a range of k between 0.02 — 0.06 h Mpc -1 reaches the 
nonlinear estimate. (Actually, the growth is somewhat 
faster in the nonlinear treatment such that the nonlinear 
estimate is reached a bit earlier). We will argue in the 
next section that the linear approximation breaks down 
for all wave numbers k < 0.06 h Mpc -1 for z < Zb- If 
one would continue the linear approximation for z = 0.8 
or z — the range of k for which the nonlinear estimate 
is reached would gradually extend towards the present 
horizon. The interpolation shown in Fig. [3] corresponds 
to a freezing of the fluctuations for k < 0.06 h Mpc -1 for 
z < z b . 



V. SCALAR POTENTIAL OF NEUTRINO 
LUMPS AND BACKREACTION 

The crucial ingredient in the analysis of neutrino fluc- 
tuations for a large neutrino-cosmon coupling f3 is the 
local scalar potential HSHil. In parallel to Newton's 
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FIG. 6. Neutrino induced gravitational potential vs. comov- 
ing wavenumber k for the varying coupling model for redshifts 
z — 0, z = 0.8, z — 1.7, z = 3.4. Green lines show the result of 
linear perturbation theory, red lines the nonlinear estimates 
for z = 0, 0.8, 1.7 for a distribution of lumps with size 15 Mpc. 



gravitational potential it determines the strength of the 
attractive force between neutrinos, being stronger than 
gravity by a factor 2/3 2 . Furthermore, the local scalar po- 
tential is directly related to the local deviation of the cos- 
mon field from its cosmological average value. The local 
cosmon field determines the local neutrino mass. Possible 
backreaction effects of an effective neutrino mass in the 
lumps being smaller than the cosmological average are di- 
rectly related to the strength of the local scalar potential. 
An investigation of the breakdown of linear perturbation 
theory may therefore concentrate on the evaluation of the 
size of non-linear effects for the values of the local cosmon 
potential and the corresponding cosmological average. 

A neutrino lump not only generates a gravitational po- 
tential in its vicinity, but also a local scalar potential 
which is associated to a local deviation S(f> of the cos- 
mon field from its cosmological average value. The force 
on neutrinos exerted by the local scalar potential is en- 
hanced by the coupling (3, cf. eq. (|16p . such that the 
quantity to be compared with the gravitational potential 
is j35(f>. As long as eq. (|17|) remains a valid approxi- 
mation, one has the simple relation S(f> = 2/3$^, cf. eq. 
(|18p . One arrives at a simple proportionality between the 
scalar potential and the neutrino induced gravitational 
potential 



= 2/3 2 $, 



(22) 



This proportionality holds both for the local potential 
5<j>i and the cosmological averaged potential 8<fi(k). Since 
the factor 2/3 2 is large, for example 2/3 2 ~ 1.5 x 10 5 for 
the model with constant /3 and m v (to) = 0.48eV, the 
scalar potential (36(f) reaches values close to one when 
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the gravitational potential is still very small. For values 
of k within the horizon one can directly infer the scalar 
potential from Fig's[3]-|n]by multiplying with a factor 
2f3 2 . We can relate the cosmological value of the scalar 
potential to the neutrino overdensity, 



P 2 A„(k)p„a 2 3/3 2 Q v (t )H 2 



A„(A) , (23) 



where the second identity holds for the present epoch 
o = l, with Hq the present Hubble parameter and il u (to) 
given by eq.(|4]). We find that /36<f> exceeds A„ for k 0.01 

Mpc -1 (model with constant j3 = —275). In this range 
strong nonlinearities in the cosmon fluctuations are ex- 
pected even before A v comes close to one. 

As we have discussed in sec lllll local values of the 
scalar potential fiSipi close to one can lead to important 
differences between a gas of neutrino lumps and a fluid 
of unclumped neutrinos. This is related to the nonlinear 
effects called "backreaction" . In presence of these nonlin- 
ear effects both the cosmological field equations for the 
background fields and the fluctuation equations for the 
modes with large wavelength are modified. One possibly 
important effect concerns the relation between the cos- 
mological average of the neutrino mass and the cosmo- 
logical average of the cosmon field. The average neutrino 
energy density, which influences H and the evolution of 
neutrino overdensities, is given by 



(24) 



where n„ is the neutrino number density. The brackets 
indicate an averaging of the local value over the horizon 
and () n is the "number- weighted average". Since m v (<j>) 
is a nonlinear function of <j>, its number weighted average 
does not equal m v {{4>)) as given by eq.©, with (0) = <j> 
the homogeneous cosmon "background field" . Typically, 
the local value of <f> inside a lump is smaller than the 
cosmological average such that the neutrino mass inside 
the lump is smaller than m v {4>). If a large part of the 
neutrinos is concentrated in lumps this leads to a "mass 
reduction factor" r m , 



m v (4>) 



< 1 



(25) 



A second effect concerns the effective cosmon-neutrino 
coupling. In eq.© (and similarly in eq.(|17p for the long 
wavelength fluctuations), one should modify the "neu- 
trino force" F v = a 2 (3p v (we neglect pressure terms) by 



F v = -a 1 



d 



(26) 



In consequence, one has to replace /3 by an averaged in- 
teraction strength (/3) as defined by 



<J3)=r p , r 



dln(r 



dhxm v {4>) 



d4> 



(27) 



In view of the nonlinearities in the field equations it is 
well conceivable that the neutrino mass inside a stabi- 
lized lump depends only very mildly on the cosmological 
average field <j>. In this case one expects a strong re- 
duction of the cosmon mediated attractive force between 
neutrino lumps, rp <C 1. This may have rather dra- 
matic consequences for the growth of fluctuations with 
long wavelength. A substantially reduced growth rate 
could effectively stop the increase of A„(fc) and $„(fc) 
for small k as soon as the smaller scale lumps form and 
most neutrinos are in lumps for which the neutrino mass 
is insensitive to the average cosmon field <fi. Whenever 
1 — r m or 1 — rp reaches values of the order one (say 
1/4) linear perturbation theory breaks down. Finally, 
also the potential derivative in eq.© has to be replaced 
by (dU/d<f>). 

The breakdown of linear perturbation theory for the 
long-wavelength modes can be seen directly by an esti- 
mate of the relative size of the nonlinear terms in the evo- 
lution equations. As an example, we consider the model 
with varying f3 and estimate the size of the nonlinear 
terms involving d/3/d(f>. For a 0— dependent j3 eq. (fl~7f is 
replaced by 



AS</> = S L + S nl 



(28) 



Sl = ~a 2 p~ 

= Sli + Sl2 



Snl — —C^Pv 



G9(0 + 50)-0(0))(i + <y„)- 



2 -dp. 



(29) 



(30) 



where we have neglected in the last line higher powers 
of 5(j>. We observe that 8<j> has the opposite sign of 8 y 
while /? and dj3/d<j) — —ft 2 are both negative. Thus in 
the linear approximation the second term Sl2 ~ S<f> par- 
tially cancels the first term Sli ~ S v . Since the common 
prefactors do not matter for a comparison of the relative 
size of the different terms, we define S — ~a 2 p v f3S. 

To evaluate the magnitude of the linear and nonlinear 
terms in Fourier space, we calculate the power spectrum 
of both terms to the first non trivial order. In Fourier 
space, the two terms Sli ( w e consider only this term for 
simplicity) and S n i can be written as 



Sli = N J 6„(k)e ikx d 3 k , 
d In /3 



(31) 



N 2 / ^(fci)<5(/>(fc 2 )e l(kl+k2)x d 3 fcid 3 /c2 



Snl = 



where N = V/(2ir) 3 is the Fourier prefactor. At any 
given wavelength fc, the two terms contribute therefore 
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as 



N-'Sl^ 8 u {k) 



We want to see for which range of redshifts and k the 
non-linear term gives a contribution similar in size to the 
linear terms. The linear approximation will clearly break 
down if this is the case. For this purpose we can take at 
a given z the linear approximation to 8 v (k) and 6<j)(k) 
and evaluate the non-linear contribution in this approx- 
imation. Only if the non-linear contribution in this ap- 
proximation remains small as compared to the linear one 
we can trust linear perturbation theory. Since by defi- 
nition of the homogeneous background the cosmological 
average of the linear fluctuations 5 v (k) and 6<j>(k) van- 
ishes we have to investigate the two-point function of the 
fluctuations of 8(f). 

The variance of the linear term is proportional to 



G L1 (k,k') = {8 v {k)8 v {k')) 



(2tt) 

V 2 



-6 3 (k + k')P ss {k) , 
(32) 

where V is the survey volume and P$$ = V\S v (k)\ 2 is the 
power spectrum. Integrating over k! we have 



(2^)3 J G Lid 3 k' = P ss (k) 



(33) 



In the Gaussian approximation one finds for the non- 
linear contribution to the two-point function to 6<f> a sim- 
ilar expression (we suppress momentarily the v subscript 
from 8 V and write momentum arguments as indices for 
clarity) 

G ni (k, k') = (^-) 2 N 2 [ (SkMkJkJA*) (34) 



x S 3 (k - ki - k 2 )S 3 (k' - fc 3 - k 4 )Y[d 3 k t . 

z=l 

Now the fourth-order term (8k 1 54>k2&k 3 54>k l ) can be de- 
composed into products of second order terms (since lin- 
ear perturbations are supposed to be Gaussian): 

(<5fci<5</>fc 2 <5fc3<5</>fc 4 ) = (Sk 1 S(t>k 2 )(8k 3 8<j)k 4 ) (35) 

+ (<5fei^fc 3 )(^fe 2 ^fe4) + {hxHk^iHk^ks) ■ 

Since {5k 1 54>k 2 ) = ^yl 8 3 {k\ + k 2 )P v ^ and similarly for 
the other pairs, we obtain, successively integrating over 
k 2 ; ^4 7 k , 



Inl = 



V 2 



(27T) 



G nl d 3 k' 



^^) 2 (2tt)- 3 I (PhlPuuikiiPMik-ki) 



+P v< f,(k 1 )P v<f) (k-k 1 )} 



(36) 



(a term proportional to S 3 (k) has been discarded since 
it vanishes for all k ^ 0). The position space variance is 
related to the power spectrum P{k) by 



a% = (8(x)8(x + r)) R = 



2tt 2 



P{k)W B 2 {k)k 2 dk (37) 



where Wn(k) is the window function of a spherical cell 
of size R in Fourier space. Eq. (|33p shows that our defi- 
nition of Il\ is related to a directly observable quantity 
and similarly for the non-linear contribution /„; . We see 
then that all factors containing the survey volume drops 
out in these quantities. (The output of CAMB gives the 
dimensionless quantities A& = 8k{k 3 V /2-k 2 ) 1 / 2 and the 
analogous expression for Stfik-) We next use the fact that 
P vu {k) and P u( j,{k) depend only on \k\ and that 



nl 



{d\n(5/d<j)) 2 
(2tt) 2 



k 1 dki / d/i 



x [Puu{k x )P^k 2 + k 2 - 2kkin) 
+ Pu^Pu^ k 2 + k 2 ~ 2kk lf i)] 



(38) 



For small k the non-linear term becomes fc-independent. 
We may now replace P vu ,P uc f, by A uv , A„0, 



P a p = 2Tr 2 A 2 a pk- 3 , 

and extract A Q ^(fc) from a CMB-code as CAMB. 
In the Newtonian approximation one has 



p*m 



k 2 
fc 4 



P V u{k) 
P vv {k) 



8tt 2 /3 
a 2 p v k 
8t: 2 [3 2 



fc 3 



-**(*) 



(39) 

(40) 
(41) 



We can therefore express I n i in terms of the cosmological 
neutrino induced Newtonian potential <!>„ and P vv 



I 



nl 



2(3 4 



din [3 
d4> 



dki 

IT 



dfi 



x <I>l{ki)PU^/k 2 + k 2 -2kk 1 n ) 

k 2 k 4 



_k 2 + kf- 2kkifi 
to be compared with Ii,\{k) 



(k 2 



k\ 



2kk X [if 



(42) 



P vv {k). For the grow- 
ing coupling with (9 In j3 / 'd<f)) 2 = f3 2 the nonlinear term 
typically reaches the same order as the linear term if <S> U 
reaches a value j3~ 2 in a relevant momentum range. 

We plot in Fig|7] the square root of the two terms 
Ili,ItiI in arbitrary units at two epochs z = 1.1 and 
z = 1.7. (At these redshifts the quantity (38(f)(k) reaches 
first a "threshold value" (38<p = 0.01 or /38(f) = 0.1 for 
some momentum k. This will play a role for the mod- 
ifications of the linear evolution discussed in the next 
section.) As we can see, the two terms are comparable 
already at z w 1.7, and the non linear term equals or 



12 



dominates (for k < 0.1/i/Mpc) the linear one at z ~ 1.1. 
If we include the Sl2 term (see Eq. |29|) . the linear term 
is further suppressed by the partial cancellation and the 
non-linear contribution becomes relevant even earlier. 

Further nonlinearities arise from deviations from the 
Gaussian approximation. The time evolution of the 
fluctuations induces non-vanishing three-point functions 
and connected four-point functions (bi-spectrum and tri- 
spectrum) In turn, these also influence the evolution of 
the two-point function 



(27T) 3 (5 3 (fc 



(43) 



k') a 4 p 2 /3 4 



V- 



ihl + IL2 + Inl + ■■■) 



and add further nonlinearities to We conclude that 
the cosmological fluctuations are already highly non- 
linear at z = 1.1. Nonlinearities presumably start to play 
an important role somewhere in the range z = 1.5 — 2. 
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FIG. 7. Comparison of the linear (full lines) and nonlinear 
(dashed lines) terms at two epochs, zi = 1.1 (thick 

colored lines) and Zb = 1.7 (thin gray lines). The dotted 
thin line is the combination of the two linear terms LI and 
L2. For z = 1.1 the linear approximation breaks down for all 
k < 0.2ft/Mpc. 



VI. MATCHING THE LINEAR AND 
NONLINEAR EVOLUTION 

Despite the present uncertainties in the estimate of the 
neutrino induced gravitational potential, one would like 
to get a first idea about the possible impact of neutrino 
lumps on the CMB. In a sense, we would like to esti- 
mate the size of the neutrino lump induced nonlinear 
Rees-Sciama effect. In order to implement the nonlinear 
features of the evolution of neutrino fluctuations in a sim- 
ple way in CMB-codes as CAMB or CMBEASY, we will 
work with simple bounds on the fluctuations that should 



reflect the dominant features of the nonlinear effects. Of 
course, this can yield at best an idea about a possible 
size of the effect. For every redshift z these bounds have 
to match the linear evolution in its range of validity with 
the information and bounds available for wave lengths in 
the nonlinear range. 

In order to employ some reasonable procedure of 
matching the linear and the nonlinear gravitational po- 
tentials, it is useful to consider the time evolution in 
terms of the physical wavenumber q. The latter is re- 
lated to the comoving scale k via k = aq. In the linear 
regime for z > z n \ we obtain for the density contribution 
of the gravitational potential: 



(44) 



where A„(q, z) changes according to the linear equations. 

We will assume that at a characteristic redshift z n i the 
non-linear effects become relevant for cosmological aver- 
ages and stop the further growth of the neutrino induced 
gravitational potential except for a change in the neu- 
trino mass, 



ml (z) 



®v{q,z n i) 



(45) 



m v {z n i) 

This corresponds to a situation where the further infall 
of neutrinos into lumps of size ix jq has a negligible ef- 
fect on the gravitational potential. We take for m v (z n i) 
the cosmological neutrino mass, while m£P(z) denotes 
the average neutrino mass in the lumps. For comoving 
momenta the gravitational potential (|45j) becomes 

'a n i\ 3 / 2 m2\z) 



(k- 



(46) 



m v {z n i) 

where the factor a -3 in <I>„ 2 (fc) is due to the defini- 
tions of the green functions (it translates the 6— function 
from q-space to k-space or, equivalently, accounts for 
the transition from physical to comoving volume). For 
$>u(q) ~ q~ a this results in a time evolution 



$ u (k,z)~m^ l \z)a a - 3 / 2 



(47) 



We will discuss two different scenarios. For the first 
one (bound 1) the backreaction effects freeze the neutrino 
mass inside the lumps, mi^(z) ~ m v {z n {). For a = 2 this 
results in 



<f> v (k,z) 



I + z n i p„(z n i)A„(z n i(k)) 
l + z 2fc 2 (l + z nJ ) 2 



(48) 



while for a — 3 /2 (cf. cq |2T)| the first factor is can- 
celled and $„(£;, z) becomes independent of z. The 
second scenario (bound 2) takes for m v "\z) the cos- 
mological neutrino mass, which multiplies the es- 
timate of the first scenario for $„ by an oscil- 
latory factor m v (z) / m v (z n i) . The numerical im- 
plementation of bound 1 corresponds to A(fc, z) — 



13 



(a b /a) 7 / 2 p l/ (z b )/p l/ (z)A(k,z b ); bound2 is given by 
A(k,z) = (a b /a) 7/2 p u (z b )/ p v (z)A(k,z b )m u (z)/m lJ (z b ). 
The truth is probably somewhere in between. 

At this stage we can devise simple strategies to imple- 
ment the non-linear effects in a CMB-code. Of course, 
this can only be used for order of magnitude estimates 
of the effects of neutrino lumps, not for the computa- 
tion of a detailed spectrum. Our general approach is 
to follow the linear evolution until its breakdown at red- 
shift z n i . Subsequently we switch to the non-linear evolu- 
tion according to eq. (|4"5j) . employing either the scenario 

with constant rriu (z) — rn u {z n i) or with varying neu- 
trino mass m v (z) according to the cosmological solution. 
For the non-linear evolution we put the velocity contri- 
bution to $„ to zero. In general, one may have different 
z n i for different k~ modes. We therefore have to device a 
reasonable way to determine z n \. 

Consider first the case that the dominant non-linearity 
arises from the neutrino overdensity. In this case we 
choose a characteristic A c J %t such that the linear theory 
breaks down for A v (k) > A£"'. We determine z n i(k) 
by the condition A l/ (z n i(k), k) — A c J %t . Thus each mode 
is switched separately to the non-linear evolution when 

A v (k) reaches the critical value. In practice we employ 
A cnt _ jy 4 

As we have discussed in sect|Vl neutrino overdensities 
close to one are not the only source of nonlinearity. We 
should incorporate the large backreaction effect discussed 
in the preceding section. Its dominant role is to reduce 
the effective interaction ((3) and the effective neutrino 
mass {m v ) n "seen" by the long-distance modes. Both ef- 
fects reduce the growth of fluctuations as compared to the 
linear approximation. One possibility would be the use of 
"modified linear equations" for which m„((/>) and P{4>) are 
replaced by effective fc-dependent functions m„ )e ff(<^, k), 
/3 c ff((j>, k). An estimate of the effective functions requires 
an understanding of the response of neutrino lumps of 
size na/k to variations of the "outside cosmon field" (f>, 
i.e. dm Vye ff(4>,k)/d(j). We assume here that this response 
is much smaller than for the case of an unclumped neu- 
trino fluid. In consequence, the reduction of the growth 
rate for neutrino fluctuations with small k is substantial, 
which justifies the use of eq.(|45|) for z < z n i. We take here 
the crude approximation that the backreaction effects set 
in at some redshift z b for all modes with k < k b . Thus for 
all modes with k < k b we take z n i = z b independently of 
k. (This holds unless A c J %t has been reached even earlier 
for some fc-modes. For these modes the "A£ ri *-trigger" 
results in z n i(k) > z b .). 

In practice we implement the backreaction effect in the 
CMB code by a "backreaction bound" for the growth 
of fluctuations, which acts in addition to the bound 
A„(fc) < A cr ; t . For this purpose we select a scale k b 
and follow the evolution of (35(f>(k b ) until a critical value 
?7 cr it is reached. At this value of z = z b we switch to 
the non- linear evolution for all modes with k < k b . For 



the modes with k > k b the backreaction bound is not 
effective such that their growth is limited only by the 
condition A„(fc) < A cr j t . For k < k b , however, the end 
of the linear evolution and the approximate freezing of 
the growth of neutrino fluctuations according to eq.(05]) 
is now set by the nonlinearity at the scale k b , and no 
longer individually for every scale k. This reflects the 
overall picture that backreaction plays only a minor role 
for k > k b , while it changes the behavior drastically for 
k < k b . In other words, our procedure mirrors an effec- 
tive fc-dependent coupling /3(fe, z) which is substantially 
smaller than f3 for z < z b and k < k b . We note that 
a strong backreaction effect would also affect the back- 
ground evolution for z < z b by a modified effective /?. 
We have not included such an effect yet. As long as 
the evolution of <j> remains very slow also in presence of 
the backreaction, the overall cosmology for z < z b , for 
which dark energy behaves very close to a cosmological 
constant, will not be modified much. 

For an estimate of an appropriate choice of k b we plot 
in Fig|8]the redshift z b (k) for which the linear evolution of 
j38(f)(k), constrained by the bound A v {k) < A cr it = 1/4, 
reaches a given value /3S(f>(k) — 0.001,0.01 or 0.1. We 



2.5 




k [h/Mpc] 

FIG. 8. Redshift z b at which the linear evolution of flS<f> 
reaches the limit value 10 -3 , 10 -2 , 10 _1 respectively, for the 
varying coupling model. 



observe that a whole range in k reaches the given value 
of (38(j)(k) almost simultaneously. This can be under- 
stood by a glance at FiglHl where for a whole range of 
k the gravitational potential $> v (k) reaches a "critical 
slope" $„, cr it/fc 2 at z = 1.7. (Since 8<f>(k) = 2£ 2 $„(fc) 
we can use Fig J6] for the time evolution of 5<p(k), pro- 
vided one rescales for every z the units on the vertical 
axes by a factor 2/3 2 (z).) We choose k b at the right 
end of the plateau in FigJH whereupon the precise choice 
within the plateau does not matter. In the varying cou- 
pling model we choose z b = 1.9(1.7,1.1) and k b — 5 x 
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1(T 2 (2.3 x 10~ 2 , 5 x 1(T 3 ) for r) crit = 10~ 3 (10~ 2 , 10" 1 ) 
respectively. In the constant coupling model Zb = 1.4 (1) 
and k b = 3x 1(T 2 (3.3 x 1CT 3 ) for r] crU = 10~ 2 (lO" 1 ) 
respectively. For the choice of the value j35(f)(k) = r] CT n, 
for which the backreaction stops the growth of fluctu- 
ations with k < kb, we recall that the cosmological 
value of (3Scj)(k) is related to the local value of the scalar 
potential of individual lumps with radius Rf — I by 
(3S4>(k) — y c p5(f)W (I — na/k) by the dilution factor 7 C in 
the range 10 -3 — 10~ 2 . Since the neutrino mass within 
the lumps differs strongly from the cosmological neutrino 
mass for |/3<5(/>W(Z)| = 1, a typical value of 77 cr it should be 
in the range of 7 C . We will display results for two choices, 
?7crit = 0.01 and r] CI it = 0.001. Furthermore, the inter- 
polation becomes very smooth if at the redshift where 
(3S(j){kb) reaches r] CTit also A v (kb) reaches A crit . This 
is realized for the pairs (A crit = 1/4, 77 crit = 0.01) and 
(A cr it = 10 , rjcrit = 10~ 3 ), that we will use for our es- 
timate of the CMB anisotropics for the growing coupling 
model. 

A typical evolution of fluctuations according to our 
switch to the non-linear evolution is shown in Figj9l 
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FIG. 9. Evolution of perturbations from the bounded linear 
code, for a scale k = 0.005 1/Mpc. The bounds are set by 
Acrit = l/4,77crit = 0.01, k b = 0.023, 1 + 2 6 = 2.7. For this 
k scale 7y cr it is reached before A cr it- Bound 1 and bound 2 
correspond to the conditions given above. 



VII. IMPACT OF NEUTRINO LUMPS ON CMB 

The impact of the neutrino spectrum on the CMB- 
anisotropies strongly depends on the evolution of the fluc- 
tuations in the nonlinear range. We demonstrate this for 
the growing coupling model with a present average neu- 
trino mass m v {to) = 0.46eV in FigJTUJ This figure shows 
the spectrum of anisotropics for various approximations: 



1. Linear perturbation theory (blue dot-dashed line) 
shows drastic deviations from the observed spec- 
trum for I < 200. This result agrees with [l5|. As 
we have argued, linear perturbation theory is not 
appropriate, however. 

2. Imposing the bound A J/ (/c) < 0.1, together with 
the backreaction bound (35(f>(kb) < 10~ 3 for fc 6 = 
5 x 10~ 2 and using the bound 1 for the evolution 
at z < z n i, we find the green dashed line. The 
neutrino fluctuations leave almost no imprint on 
the spectrum. Deviations from the ACDM model 
result primarily from the presence of early dark en- 
ergy. 

3. For a less restrictive bound A„(fc) < 0.25, 
(35(f>(kb) < 10 -2 the neutrino induced gravitational 
potential is somewhat larger (cf. the upper vs the 
lower dotted blue line (bound 1) in Fig|3j)) . Now 
a signal of the neutrino fluctuations is visible as 
an enhancement of the anisotropics in the range 
10 < I < 80, as shown in the red solid line in FigflUl 
This gives an indication of the size of the ISW - or 
nonlinear Rees-Sciama effect for a neutrino induced 
gravitational potential of the size shown in Fig|3l 
The smoothness of this effect is a result of the par- 
ticular type of bounds employed. Other implemen- 
tations of the nonlinear neutrino fluctuations can 
lead to an oscillatory behavior. 

4. A much stronger ISW effect would result for an 
increased neutrino induced gravitational potential. 
The purple dotted line shows the spectrum for the 
bound f38<j){kb) < 0.1. Strong discrepancies with 
observations would occur for / < 80. In view of 
the estimate of the backreaction effect in sect I VII 
this curve cannot be trusted, however. The em- 
ployed linear approximation is invalid within the 
range of redshifts during which the large fluctua- 
tions responsible for the ISW-effect grow. 

For all curves we use the cosmological parameters 
shown in the first column of table (JXJ) . These parame- 



growing coupling 


fixed coupling 


m„(to) = 0.46 eV 


m v (t ) = 0.48 eV 


a = 10 
4>t = -0.005 


a = 10 
P = -275 


O^o = 0.67 
Q c0 = 0.26 
fi„o = 2.82 • 10~ 2 
ftbo = 4.14 ■ 10~ 2 
tfo = 72.9 


= 0.65 
n c0 = 0.276 
O,„o = 3.12 • 10~ 2 
Q b0 = 4.34 ■ 10~ 2 
Ho = 71.3 



TABLE I. Cosmological parameters for a varying P{4>) model 
and a constant P model. 
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ters are modified as compared to the ones for the ACDM 
model due to the presence of early dark energy. For 
larger values of a, resulting in a reduced early dark en- 
ergy, they would be closer to the ACDM model. We have 
not attempted an optimization of parameters in this pa- 
per. Deviations of the shown spectrum from the ACDM 
model for I > 200 can be reduced by an improved pa- 
rameter choice. Restrictions on models arising from the 
CMB spectrum for I > 200 or the corresponding power 
spectrum of matter fluctuations for k > 3 x 10~ 2 Mpc -1 
mainly concern a bound on the parameter a and are not 
addressed in this paper. 

The CMB-anisotropies are sensitive to the neutrino 
mass and the form of the neutrino-cosmon coupling. In 
FiglTT]we plot results for similar bounds as in Fig[TUl but 
now for a fixed coupling f3 — —275. The cosmological pa- 
rameters used for this figure are displayed in the second 
column of Tab|T|. 
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FIG. 11. Angular power spectrum of CMB anisotropies for 
the constant coupling model, /3 = —275. 
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FIG. 10. Angular power spectrum of CMB anisotropies for 
the varying coupling model. The lines correspond to differ- 
ent approximations for the nonlinear neutrino fluctuations as 
explained in the text. 



VIII. CONCLUSIONS 

We have attempted to estimate the size of the ISW- 
effect for the CMB spectrum which results from the for- 
mation of large scale neutrino lumps in models of growing 
neutrino quintessence. Many of our general considera- 
tions also apply for other models for which the coupling 
between neutrinos and a light scalar field substantially 
exceeds gravitational strength, as for the mass varying 
neutrino models of 0-@|- Our most central conclusion 
seems quite robust: linear perturbation theory cannot be 
trusted for an estimate of the CMB spectrum for these 



models. This often concerns already the onset of observ- 
able differences from the ACDM model, since the neu- 
trino fluctuations which dominate the ISW-effect in the 
linear approximation actually correspond to a range of 
wave numbers k for which the linear approximation is no 
longer valid. 

For an estimate of the CMB-spectrum for I < 100 
in growing neutrino quintessence one needs to solve a 
rather complex nonlinear problem, similar to the esti- 
mate of the nonlinear Rees-Sciama effect for I > 1500 
arising from dark matter lumps. In contrast to the dark 
matter lumps, N-body simulations or reliable analytic re- 
sults for the neutrino-lump-induced gravitational poten- 
tial $> v (k) in the relevant range k < 0.04 h Mpc -1 are not 
yet available. The challenge for such methods consists in 
an appropriate treatment of a locally varying scalar field 
in the presence of possibly important backreaction ef- 
fects. Waiting for more accurate results for $„(fc) and 
its evolution with redshift, we attempt an estimate of 
the plausible overall size of the ISW effect by imposing 
reasonable bounds on the growth of fluctuations. 

For the varying coupling model we learn from Fig llOl 
that an observable ISW-effect for Z < 100 seems plausi- 
ble, perhaps dominating around ! w 40. It seems well 
conceivable that this effect includes oscillatory features, 
reflecting oscillations in the cosmological neutrino energy 
density as visible in Fig|2j One may even speculate that 
this could be responsible for small "enhancements" and 
"depressions" in the observed spectrum for this range 
of I. The size and details of the neutrino- induced ISW- 
effect depend on the particular model and its parameters. 
Unfortunately, a detailed parameter estimate seems pre- 
mature in view of the large uncertainties affecting the av- 
erage gravitational potential of the neutrino lumps. Nev- 
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ertheless, it may be possible to exclude certain ranges of 
models and parameters if it is plausible that the devia- 
tions from the linear behavior in the relevant ranges of k 
and z are under control. 

A general lesson concerns the presence of strong back- 
reaction effects within growing neutrino quintessence. A 
gas of neutrino lumps behaves differently from an ap- 
proximately homogeneous fluid of neutrinos. The mi- 
croscopic values m v (4>) and (3{4>) of the neutrino mass 
and cosmon-neutrino coupling get replaced by an effec- 
tive mass and coupling once the neutrino lumps have 
become sufficiently nonlinear. This seems to be one of 
the first examples of a realistic large backreaction ef- 
fect which modifies the parameters of the homogeneous 
and isotropic background field equations by order one 
effects. A self consistent treatment beyond linear pertur- 
bation theory will have to take them into account. As a 
whole, the findings of this paper are encouraging for more 
detailed investigations of growing neutrino quintessence. 
It seems that such models could be consistent with ob- 
servations for appropriate ranges of parameters. Once 
the true evolution of the neutrino-induced cosmological 
gravitational potential z) is better constrained, in- 

teresting model tests become feasible through neutrino- 
induced modifications of the CMB-spectrum for I < 100, 
an enhanced correlation between matter and photon fluc- 
tuations or an enhanced large scale bulk flow. The pos- 
sible existence of large scale neutrino lumps could offer 
new surprises for cosmology. 

We thank David F. Mota for providing us with a first 
version of CAMB code including mass varying neutrinos. 
The work of NW is supported by the Humboldt Founda- 
tion. 
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APPENDIX 

LINEAR PERTURBATIONS WITHIN GENERAL 
GROWING NEUTRINO MODELS 

We recall here the linear perturbation equations for 
growing neutrino quintessence, following the work pre- 
sented in [12J, to which we refer for further details. The 
evolution equations for linear perturbations (in Fourier 
space), in Newtonian gauge (in which the non diagonal 
metric perturbations are fixed to zero) 37J, read for the 
growing neutrino scenario: 



= 3"H(w0 - c|)(5 

- [(1 - 3tu tf )(5* - (1 - 3cl)5„] 
P<t> 

- (l + u^)(fctv + 3$') 

+ EJL {1 _ 3^) ( m sct>' + 6^ 



, (49) 



- (1 + u>„)(fct;„ + 3$') - - 3«v )£</>' 



(50) 



The equations for the density contrasts 5i(k) — 
y J Si(x)exp(— ik ■ x)d 3 x (defined as the Fourier trans- 
formation of the local density perturbation <5j(x) = 
8pi{x)l Pi{x) over a volume V) involve the velocity per- 
turbations, which evolve according to 

i>0 = —H(l - iw^Vtj, - /3(4>)<fi'(l - 3w v )—v<f, 



1 + W<h V 1 + Wtf, 



W<j, 

2 



3 1 + w<f, 



i+^0 



^ = {l-ZwJtfWtf -U)v v - 



1 + w v 



kct 



1 T 2 , ^ 

fcv — — k- T^T v 



1 + w v 3 1 + w v 

1 — 3u)^ 



- WW 1 

1 + w„ 

As usual, the gravitational potentials obey 

„2 



2k 2 M 2 



w a p a n Ta 



(52) 



(53) 



(54) 



where 7Tt q is the anisotropic stress for the species a and 
the sound velocities are defined by c 2 = Spi/Spi. The 
perturbed pressure for <f> is 

8 V4> = ^60' - -?U' 2 - (55) 

and the anisotropic stress 7Tj\ = 0, as in uncoupled 
quintessence, since the coupling is treated as an external 
source in the Einstein equations. The linear perturbation 
of the cosmon, 5(f), is related to via 

8<j> = (fi'v^/k , 

+a 2 /3((/))(p L , - 3p„)] V4, . 

Note that 8<f> can equivalently be obtained as the solution 
of the perturbed Klein Gordon equation: 
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+ 2a 2 ^-^ = -a 2 \-P((t))p v 8 u (l-3cl) (57) 
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The evolution of neutrinos requires solving the Boltz- 
mann equation in the case in which an interaction be- 
tween neutrinos and the cosmon is present 38] . The first 
order Boltzmann equation written in Newtonian gauge 
reads [H 



or e 



din/, 



= i-(k • n)fc 



ps ' dlnq 
a 2 m 2 <91nm„ din fo 
q 2 d(f> din g 



$' - i-(k ■ n)* 



f(x\T,q,n j ) = f (q) [l + ^ ps (x\ r, q, nj)] 



(58) 



(59) 



*5> and $ are the metric perturbations, x l the spatial co- 
ordinates, t is the conformal time, q = ap = qh is the 
comoving 3-momentum, e = e(0) = \J q 2 + m u (4>) 2 a 2 , f 
is the phase space distribution and fo its zeroth-order 
term (Fermi-Dirac distribution). 

The Boltzmann hierarchy for neutrinos, obtained ex- 
panding the perturbation ty ps in a Legendre series can 
be written in Newtonian gauge as 



1 qk , T „ _ . ek _ d In fn 



(60) 



Vj/' , = - 

(2Z + l)e 



where 38, 3 
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This allows us to calculate the perturbed energy and pres- 
sure as well as the shear for neutrinos: 
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The anisotropic stress is related to the shear via 7Tt„ 
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~T~{Pv + Vv) G arL d in our case 
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Note also that the unperturbed neutrino density and 
pressure read 
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